NASA/TM— 2007-214927 



Progressive Fracture of Laminated Composite 


Pascalis K. Gotsis 

Technological Educational Institute of Serres, Terma Magnisias, Serres 
Christos C. Chamis 

Glenn Research Center, Cleveland, Ohio 
Kostantinos David 

Technological Educational Institute of Serres, Terma Magnisias, Serres 
Frank Abdi 

Alpha STAR Corporation, Long Beach, California 


Stiffened 



August 2007 


NASA STI Program ... in Profile 


Since its founding, NASA has been dedicated to the 
advancement of aeronautics and space science. The 
NASA Scientific and Technical Information (STI) 
program plays a key part in helping NASA maintain 
this important role. 

The NASA STI Program operates under the auspices 
of the Agency Chief Information Officer. It collects, 
organizes, provides for archiving, and disseminates 
NASA’s STI. The NASA STI program provides access 
to the NASA Aeronautics and Space Database and its 
public interface, the NASA Technical Reports Server, 
thus providing one of the largest collections of 
aeronautical and space science STI in the world. 
Results are published in both non-NASA channels and 
by NASA in the NASA STI Report Series, which 
includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant phase 
of research that present the results of NASA 
programs and include extensive data or theoretical 
analysis. Includes compilations of significant 
scientific and technical data and information 
deemed to be of continuing reference value. 
NASA counterpart of peer-reviewed formal 
professional papers but has less stringent 
limitations on manuscript length and extent of 
graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies that 
contain minimal annotation. Does not contain 
extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, often 
concerned with subjects having substantial 
public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific and 
technical material pertinent to NASA’s mission. 

Specialized services also include creating custom 

thesauri, building customized databases, organizing 

and publishing research results. 

For more information about the NASA STI 

program, see the following: 

• Access the NASA STI program home page at 
http://www.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA STI Help Desk 
at 301-621-0134 

• Telephone the NASA STI Help Desk at 
301-621-0390 

• Write to: 

NASA Center for AeroSpace Information (CASI) 
7115 Standard Drive 
Hanover, MD 21076-1320 


NASA/TM— 2007-214927 



Progressive Fracture of Laminated Composite 


PascalisK. Gotsis 

Technological Educational Institute of Serres, Terma Magnisias, Serres 
Christos C. Chamis 

Glenn Research Center, Cleveland, Ohio 
Kostantinos David 

Technological Educational Institute of Serres, Terma Magnisias, Serres 
Frank Abdi 

Alpha STAR Corporation, Long Beach, California 


National Aeronautics and 
Space Administration 

Glenn Research Center 
Cleveland, Ohio 44135 


Stiffened 



August 2007 


Acknowledgments 


The first author would like to thank the Research Committee of Technological Education Institute of Serres, Greece, for 

funding this research project. 


This work was sponsored by the Fundamental Aeronautics Program 
at the NASA Glenn Research Center. 


Level of Review. This material has been technically reviewed by technical management. 


Available from 


NASA Center for Aerospace Information 
7115 Standard Drive 
Hanover, MD 21076-1320 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22161 


Available electronically at http://gltrs.grc.nasa.gov 


Progressive Fracture of Laminated Composite Stiffened Plate 

Pascalis K. Gotsis 

Technological Educational Institute of Serres 
Terma Magnisias, Serres 62120 

Christos C. Chamis 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Kostantinos David 

Technological Educational Institute of Serres 
Terma Magnisias, Serres 62120 

Frank Abdi 

Alpha STAR Corporation 
Long Beach, California 90804 

Abstract 

Laminated fiber-reinforced composite stiffened plate with [0/90/±45]s plies made of S-Glass/epoxy 
are evaluated via computational simulation to study damage and fracture progression. The loads are 
pressure and temperature which varies from 21 to 65.5 °C (case I) and from 143.3 to 21 °C (case II). An 
integrated computer code is used for the simulation of the damage progression. Results show that damage 
initiation begins at low load level, with matrix cracking at the 0° (bottom and top) plies, fiber fracture at 
the bottom (0°) ply and interply delamination at the top (0°) ply. Increasing the applied pressure, the 
damage growth is expended resulting in fracture through the thickness of the structure. At this stage, 

90 percent of the plies damage at applied pressure 15.306 MPa for the case I and 15.036 MPa for the 
case II. After this stage the cracks propagate rapidly and the structure collapses. 

1. Introduction 

The aircraft, marine and automotive industries use stiffened composite plates because of their low 
weight, high stiffness, and stability. Design considerations with regard to the durability stiffened plates 
require a priori evaluation of the damage initiation and propagation mechanisms under expected service 
loads. Concerns of the safety and survivability of critical components require quantifications of the 
composite structural damage tolerance during overloads. Characteristic flexibilities in the tailoring of 
composite structures make them more versatile for fulfilling structural design requirements. However, 
these same design flexibilities render the assessment of composite structural response and durability more 
complex, prolonging the design and certification process and adding to the cost of the final product. It is 
difficult to evaluate composite structures because of the complexities in predicting their overall congruity 
and performance, especially when structural degradation and damage propagation occur. The prediction 
of damage initiation, damage growth, and propagation to fracture are important in evaluating the load- 
carrying capacity, damage tolerance, safety, and reliability of composite structures. Quantification of the 
structural fracture resistance is also required to evaluate the durability /life of composite structures. The 
most effective way to obtain this quantification is through integrated computer codes that couple 
composite mechanics with structural analysis and damage progression models. GENOA computer code 
used in this study. The code is an updated and improved version of CODSTRAN (Chamis et al., 1996). 

The simulation of progressive fracture has been verified to be in reasonable agreement with 
experimental data from tensile coupon test on graphite/epoxy laminates (Irvin and Ginty, 1986) damage 
progression in carbon fiber reinforced plastic I-beams (Hang and Minnetyan, 1998) and carbon/carbon 
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composite plate specimens subjected to three-point bending (Solinsky et al., 2007). A variety of laminated 
fiber-reinforced composite structures used to simulate the damage progression and fracture, such as: 
stiffened adhesively bonded composite structures (Gotsis et al., 1995), damage progression in bolted 
composite structures (Chamis et al., 1995) and damage tolerance of composite pressurized thin shell 
structures (Gotsis et al., 1996). 

The purpose of this paper is to perform computational simulation to S-Glass/epoxy (0/90/±45)s 
laminated reinforced-fibers composite stiffened plate, subjected to thermo-mechanical loads, in order to 
predict the damage progression, fracture thought the thickness and propagation to final fracture (or 
collapse) of the structure. 


2. Methodology 

The computational simulation is performed by coupling three modules: (1) composite mechanics, 

(2) finite element analysis, and (3) damage progression tracking. The damage progression module relies 
on the composite mechanics code ICAN (Murthy and Chamis, 1986) for composite micromechanics, 
macromechanics and laminate analysis, and calls a finite element analysis module that uses anisotropic 
thick shell elements to model laminated fiber-reinforced composite structures. The finite element module 
is based on the mixed finite element method (Nakazawa et al., 1987). By supplying the boundary 
conditions, the type of analysis desired, the applied loads, and the laminate properties, the module 
performs the structural analysis. In addition the finite element module provides the computed stress 
resultants to composite mechanics module, which continuously computes the developed ply stresses for 
each ply and checks for ply failure. 

A computational simulation cycle begins with the definition of constituent properties from a material 
databa nk . Composite ply properties are computed by the composite mechanics module. The composite 
mechanics module also computes through-the-thickness structural properties of each laminate. The finite 
element analysis module accepts the composite properties that are computed by the composite mechanics 
module at each node and performs the analysis for a load increment. After an incremental finite element 
analysis, the computed generalized nodal force resultants and deformations are supplied to the composite 
mechanics module that evaluates the nature and amount of local damage, if any, in the plies of the 
composite laminate. Individual ply failure modes are determined using failure criteria associated with a) 
the negative and positive limits of the six ply stresses components (the in plane stresses (a m , a/12. a/22) and 
the interlaminar stresses (a/33, a/13. a/23)), b) a modified distortion energy (MDE) combined stress failure 
criterion, and c) interply delamination due to relative rotation of the plies. 

The generalized stress-strain relationship for each node is revised according to the composite damage 
evaluated by the composite mechanics module after each finite element analysis. The model is 
automatically updated with a new finite element mesh and properties, and the structure is reanalyzed for 
further deformation and damage. If ply failure criteria indicate new or additional damage during a load 
increment, the damage tracking module degrades the composite properties affected by the damage and 
reanalyzes the structure under the same load. When there is no indication of further damage under a load, 
the structure is considered to be in equilibrium. Subsequently, another load increment is applied leading 
to possible damage growth, accumulation, or propagation. In the computational simulation cases 
presented in this paper, analysis is stopped when commencement of the damage propagation phase is 
indicated by laminate fracture. Laminate fracture is predicted when major principal failure criteria are met 
for all plies at a node and for all nodes at the structural fracture path until the structure separated into two 
or more parts. 


3. Laminated Fiber-Reinforced Composite Structure 

The structure used for this investigation is a continuous fiber-reinforced laminated composite 
stiffened plate (fig. 1). The geometry of the structure has length equal to 254 mm (10 in.) and width equal 
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Figure 1 . — A laminated fiber-reinforced 
composite stiffened plate. 


to 127 mm (5 in.). The panel and the stiffened bands are made of the same polymer matrix composite 
materials S-Glass fiber and high strength (IMHS) epoxy matrix. The S-Glass fiber and epoxy matrix 
properties are given in tables 1 and 2. The fiber volume ratio is 0.60. The laminate configuration is 
[0/90/±45] s . The skin laminate of the panel consists of eight 0.3175 mm (0.0125 in.) plies, resulting in a 
composite thickness of 2.54 mm (0.10 in.). The stiffened bands well bonded to the bottom surface of the 
panel, figure 1. The stiffened band has laminate configuration [0/90/±45] s , with ply thickness equal to 
0.1587 mm (0.00625 in.) resulting in a total thickness equal to 1.27 mm (0.05 in.). The 0° plies are in the 
X-axis, figure 1. Regarding the ply numbering convention the first ply lays at the bottom surface of the 
plate, while the last ply lies at the top surface of the plate. The loads are the pressure normal to the plate, 
and a uniform temperature, which is increased from the room temperature 21 to 65.5 °C and remains 
constant. For comparison reasons another case was examined with the same pressure and the temperature 
to change from 143.3 °C to room temperature (21 °C). In the second case residual stresses are taking into 
account due to the heating and cooling of the fiber composites due to the fabrication process. The 
boundary conditions are fixed supported. 


TABLE 1— S-GLASS FIBER PROPERTIES 


Number of fibers per end 

Fiber diameter, mm (in.) 

Nonnal Modulus, GPa (psi) 

Longitudinal 

Transverse 

Poisson’s ratio 

Vl2 

v 23 

Shear modulus, GPa (psi) 

Gl2 

G 23 

Thermal expansion coefficient, /°C (/°F) 

Longitudinal 

Transverse 

Heat conductivity, J-m/hr/m 2 /°C (BTU-in./hr/in. 2 /°F) 

Longitudinal 

Transverse 

Heat capacity, .1 Kg C (B IT. lb T) 

Strength, MPa (ksi) 

Tensile 

Compressive 


200 

0.00762 (0.3xl0“ 3 ) 

85.5 (12.4x10") 

85.5 (12.4x10") 

0.20 

0.20 

35.67 (5.17x10") 

35.67 (5.17xl0 6 ) 

0.509xl0“ 5 (0.280xl0“ 5 ) 
0.509xl0“ 5 (0.280xl0“ 5 ) 

(5.208xl0“ 2 ) 

(5.208xl0 -2 ) 

712(0.17) 

2482 (360) 

2068.33 (300) 
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TABLE 2. — IMHS EPOXY MATRIX PROPERTIES 


Matrix density. Kg/m 3 (lb/in 3 ) 3x10' 7 (0.0443) 

Normal Modulus, GPa (ksi) 3.394 (500) 

Poisson’s ration 0.35 

Coefficient of thermal expansion, /°C (/°F) 0.7704 (0.428X10 -4 ) 

Heat conductivity, J-m/hr/m 2 /°C (BTU-in./hr/in. 2 / °F) (8.680E-03) 

Heat capacity, J/Kg/°C (BTU/lb/°F) 738 (0.250) 

Strength, MPa (ksi) 

Tensile 103.41 (15.) 

Compressive 241.3 (35.) 

Shear 89.62(13.) 

Allowable strain 

Tensile 0.02 

Compressive 0.05 

Shear 0.045 

Torsional 0.045 

Void conductivity, J-m/hr/m 2 /°C (BTU-in./hr/in. 2 /°F) 16.8 (0.225) 

Glass transition temperature, °C (°F) 216 (420) 


The finite element mesh consists of 396 elements and 442 nodes. Thick shell element was used for 
the computational simulation. 

Damage initiation, growth, accumulation, and propagation to fracture and to structural collapse were 
monitored as the pressure is gradually increased. 


4. Results and Discussion 

Discussion of the results of the laminated fiber-reinforced composite stiffened plate will follow for 
the two cases, the case I with the temperature to change from 21 to 65.5 °C, and the case II from 143.3 to 
21 °C. The pressure load is increased gradually in both cases. 

a) Damage initiation. In case I the damage started at applied pressure 0.55 MPa (80 psi). In the middle 
stiffener, matrix cracking in tension due to transverse tensile failure mode occurred at the 0° plies (the 
bottom and the top), fiber fracture due to the longitudinal compressive failure mode happened at the 
bottom ply 0° and delamination due to relative rotation at the top ply (0°). At the panel area, matrix 
cracking in tension due to transverse tensile failure mode occurred at the 0° (top) ply. 

b) In case II the damage began earlier at 0.275 MPa (40 psi). In the middle stiffener, matrix cracking 
in tension due to transverse tensile failure mode occurred at the 0° plies (the bottom and the top), fiber 
fracture due to the longitudinal compressive failure mode happened at the bottom ply 0°. At the panel 
area, matrix cracking in tension due to transverse tensile failure mode occurred at the 0° (top) ply. 

c) Damage progression, fracture and structural collapse. In general, overall structural damage may 
include individual ply damages and through-the-thickness fracture of the composite laminate. The type of 
the damage growth and the sequence of damage progression depend on the composite structure, loading, 
material properties, and hygrothermal conditions. A scalar damage variable, derived from the total 
volume of the composite material affected by the various damage mechanisms, is also evaluated as an 
indicator of the level of overall damage induced by loading. This scalar damage variable is useful for 
assessing the overall degradation of a given structure under a prescribe loading condition. The rate of 
increase in the overall damage during composite degradation may be used as a measure of the structural 
propensity for fracture. Damage is defined as the volume of the damaged plies divided by the total 
volume of the structure. Figure 2 presents a comparison of the overall damage progression histories of the 
stiffened plate for the two loading conditions, and it is plotted the percentage of the damaged plies versus the 
applied pressure. In case I fracture through the thickness occurred when the percentage of the ply damages was 
89.47 percent at applied pressure 15.306 MPa (2220 psi), followed by rapidly damage growth and additional 
fracture progression resulting in the collapse of the structure, at 25.09 MPa (3640 psi), at 90.23 percent damage. 

In case II fracture through the thickness occurred when the percentage of the ply damage was 89.79 percent at 
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Figure 2. — Damage progression under pressure and 
temperature S-Glass Epoxy (0/90/45/-45)s stiffened 
plate. 



Figure 3. — Dominant failure modes of the S-Glass 
Epoxy (0/90/45/M5)s stiffened plate under pressure 
and temperature (21 to 65.5 °C). 

applied pressure 15.03 MPa (2180 psi), figure 2. The collapse of the structure occurred at 27.57 MPa (4000 psi), 
at 90 percent damage. 

d) Ply failure modes. The plies of the laminated composite structure can be damaged by one or more 
failure modes simultaneously under the effect of the applied load. The various failure modes of the 
damaged plies of the stiffened plate, versus the applied pressure (case I) are plotted in figure 3. It is 
concluded the following: 

1) The MDE and the transverse tensile failure modes dominate the other modes. At the fracture 
initiation (P = 15.306 MPa) almost all damaged plies failed due to matrix cracking in tension, figure 3. 

2) The normal compressive and the longitudinal tensile failure modes are increased gradually and 
near the fracture initiation increase rapidly up to the collapse. At the fracture initiation, 62 percent of the 
damaged plies failed due to matrix cracking in compression and 36.9 percent failed due to fiber fracture in 
tension, figure 3. 
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Figure 4. — Fracture progression at P = 17.65 MPa 
(2560 psi). Stiffened plate S-Glass Epoxy 
(0/90/45/-45) s with temperature (21 to 65.5 °C). 


3) The relative rotation failure mode begins to increase slowly from the damage initiation 

(P = 0.55 MPa) state, and speed up close to the fracture initiation, where all damaged plies failed due to 
interply delamination, figure 3. 

4) The in-plane-shear failure mode increases gradually. At the fracture initiation stage 26 percent of the 
damaged plies failed and at the collapse of the structure 40. 1 percent failed due to the matrix cracking, 
figure 3. 

5) The transverse normal shear failure mode gradually increases up to the structural collapse. At the 
fracture initiation stage, 21.5 percent of the damaged plies failed due to the interply delamin a tion. The normal 
tensile and the transverse compressive failure modes have insignificant effect. 

An extensive damage in the structure, at P = 17.65 MPa (2560 psi), is shown in figure 4. 

5. Conclusions 

A computational simulation was used to evaluate the structural damage progression response of a S- 
Glass/epoxy [0/90/±45] s laminated fiber-reinforced composite structure underpressure and thermal load 
(case I and case II). The conclusions follows: 

1) Damage initiation began at low load level at 0.55 MPa for the case I and at 0.275 MPa for the 
case II. In both cases, at the middle stiffener, matrix cracking in tension due to transverse tensile failure 
mode occurred at the 0° ply (the bottom and top) and fiber fracture due to longitudinal compressive 
failure mode happened at the bottom (0°) ply. In the panel area matrix cracking in tension happened at the 
top ply (0°). In addition interply delamination took place at the top ply (0°) of the middle stiffener in the 
load case I. 

2) Fracture through the thickness began at 90 percent ply damages, due to matrix cracking, fiber 
breaking and interply delamination. The fracture initiation load for the case I was 15.036 MPa and for the 
case II 15.03 MPa . 

3) The fracture initiation followed by a rapidly crack growth resulting in the collapse of the structure. 
The collapse load for the case I was 25.09 MPa and for the case II 27.57 MPa. 

4) The dominant failure mode was the transverse tensile, followed by the normal compressive, 
longitudinal tensile and the relative rotation. The in-plane shear failure mode had moderate effect. 
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